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INTRODUCTION 


In 1952 a new approach to the correlation and extrapolation of creep- 
rupture data was introduced [1]. The approach, called the time-temperature 
parameter (TTP) method, combines time and temperature into a single expres- 
sion to collapse a family of isothermal curves into a single one. The re- 
sulting curve is called the ’’master curve.” The method assures that for a 
given stress level, rupture (or excessive deformation) will occur at any 
time and temperature combination which produces a value which corresponds 
to the value of the TTP at that stress. Various forms of the TTP expression 
have been introduced. Often mentioned in the literature are the Larson- 
Miller parameter introduced in 1952 [1]; the Manson-Haf erd parameter introduced 
in 1953 [2]; and the Orr-Sherby-Dorn parameter proposed in 1954 [3]. 

Another approach is the Minimum Commitment Method, suggested by Manson in 
1968 [4]. This method avoids forcing the data to fit a specified TTP. In- 
stead it allows the best fitting parameter to be used. And if none of the 
known parameters best fit the data, it is capable of developing a new param- 
eter. 

Over the years, significant consideration has been given to both physics 
and mathematics in the development of TTP’s. Special emphasis has been in 
the development of methods for the calculation of empirical constants for the 
parametric expressions, e.g., the methods presented by Conway [5]. Develop- 
ment of a functional relationship between strength and parameter values, 
i.e., an analytical representation of the master curve has also been 
pursued [6]. Simultaneously, studies have been made to assess the effective- 
ness of the parameter regarding correlation, interpolation and extrapolation 
of the experimental data (e.g., Goldhoff [7]). The present study applies 
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probabilistic design theory to creep-rupture data analysis, bringing state 
of the art methods to the design of components under creep. It focuses 
on the development of an analytical representation of the master curve, 
the calculation of the parameter constants, and the assesment of the ef- 
fectiveness of the parameters regarding correlation and extrapolation of 
the data. The Larson-Miller (LM) and the Manson-Haf erd (MH) parameters 
are used in this study, but the methods presented could apply to any of the 
tip’s. The Larson-Miller parameter is given by the expression 

P = (T + 460) (log^^t + C) (1) 

and the Manson-Haf erd parameter, 

T - T 

P = £ (2) 

i°giot - 

where T is temperature in °F, t is time in hours. C, T^, and log^^t^ 

are corresponding parameters or empirical constants established by the data. 

The goal of this study was to develop a statistical model for describing 
creep strength of a material using the TTP concept. This model should 
quantify scatter in material behavior as well as modelling error. Moreover, 
the model should be able to fit into a reliability format in which due 
consideration is given to all sources of uncertainty. Proposed herein is 
a lognormal model for creep strength and a lognormal format for the general 
reliability problem. 

CREEP STRENGTH AND THE MASTER CURVE 

Probabilistic design (or mechanical reliability) refers to the process 
of quantifying uncertainty, and then making decisions so that the risk is 
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less than an acceptable value. Risk is assumed herein to be synonymous 
with probability of failure. Failure of a component under creep can be 
defined as (a) the event that time to failure is less than the intended 
service life, or (b) the event that creep strength R is less than the applied 
stress S, i.e. , (R<S) . The latter definition is adopted here. Figure 1 
illustrates the basic problem. Both creep strength and applied stress are 
considered to be random variables. The risk or probability of failure is 
defined as 

p^ = Pr (R < S) (3) 

To evaluate p^, it is necessary to establish the statistical distribution of 
R and S. The process of development of the functional relationship between 
strength and the TTP also addresses the problem of translating the creep- 
rupture data into a statistical distribution of R. The basic relationship 
proposed here for creep strength is called the ’’exponential model,” 

Ios^qR = a + bp“ (4) 

A, B, and roi are parameters to be determined from the data. 

Note: In the case of the LM parameter the value of the parameter P used 

in the exponential model is in thousands , while for the MH parameter the 
absolute value of P is used. In the range of interest, P is a negative 
valued function which can not be raised to a rational power. 

The median curve of R given P can be established by a least squares 
analysis. The model parameters A, B and m are calculated by an iterative 
procedure. Note that if m is considered a constant, best linear unbiased 
estimates of the parameters A and B can be calculated by applying a least 



Creep Rupture Strength 
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Figure 1. An Illustration of the Basic Creep Reliability Problem. 
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squares analysis to a linearized form of the exponential model [8]. Set 

Y = log^gR (5) 

X = P“ (.6) 

then the exponential model becomes 

Y = A + BX (7) 

The creep-rupture data (R^, P^) can be expressed as (y^, x^) , i = l,n 
using Eqs. 5 and 6. The least squares line is defined here as those values 
of A, B and m (denoted as A, B, and m) which minimize the sample standard 
deviation s. The sample standard deviation is 



where 

t h 

y. = value of i point with coordinates (y^, x^) , 

Y. = A + Bx. 

1 1 

n - number of data points. 

ESTIMATION OE THE EXPONENTIAL MODEL PARAMETERS AND EVALUATION OP THE TTP 
CONSTANTS 

Evaluation of the TTP constants (e.g., C, T , log ^t ) and the exponen- 

a lu a 

tial model parameters (A, B, and m) is achieved using an iterative process. 

The basic criterion is that the best estimates of all of these material 
dependent constants are those which minimize the sample standard deviation, s. 
The numerical procedure is as follows (using Larson-Miller as an example) : 
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1. Choose an initial value of C. 

2* Choose an initial value of m, 

3. Use the simple least squares analysis to compute A and 

4. Compute s (Eq. 8). 

5. Repeat steps 2 through 4 until an m is found (corresponding to your 
initial choice of C) which minimizes s. 

6. Going back to step 1, repeat the process using another value of C. 
Finally, this procedure produces a minimum s as a function of C. 

7. The "best" estimate of C is defined as that value which corresponds 

A ^ A 

to the minimum value of s. The corresponding values m. A, and B 
are estimates of m. A, and B respectively. 

This procedure is easily extended to the case of two (or more) constants. 

For the MH parameter, steps 1, 6, and 7 are extended to accomodate two 
constants. 

Finally, the "optimW* master curve is 

log^^R = A + BP”" (9) 

Examples of this analysis applied to a sample of n = 95 Incoloy 625 data 

are provided in Figs. 2 (for LM) and 3 (for MH) . The coherence of the data 

to the master curve is measured by C , the coefficient of variation of R 

K 

(proportional to s' as given by Eq. 11 below) which defines the vertical 
scatter. The fact that the MH analysis produces a slightly lower suggests 
a slightly better fit, . . . not surprising in view of the additional con- 
stant in the MH parameter. 

Table 1 is a summary of the results obtained from an exponential model 
analysis for five materials. The significance of these results lies in 



Fig. 2. Example of Proposed Analysis Procedure Applied to Incoloy 625 Creep-Rupture Data Using 
Larson-MIller Parameter. 
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Fig. 3. Example of Proposed Analysis Procedure Applied to Incoloy 625 Creep-Rupture Data Using 
the Manson-Haferd Parameter. 
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TABLE 1. Sunimary of Exponential Model Analysis.* 
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the demonstration that typical values of are less than 10%. This value 
includes both scatter in material behavior and modelling error associated 
with the use of a TTP. Note that in comparison, scatter in room temperature 
static properties for most metallic materials is characterized by coefficients 
of variation of 5 to 10%. Finally, note that the R-P data for all the 
materials of Table 1, when plotted, would exhibit the same coherence to 
the master curve as those of Figs. 2 and 3. 

A STATISTICAL MODEL FOR CREEP STRENGTH 

Having estimated the model parameters we proceed to construct a statis- 
tical model of creep strength R. In general, for probabilistic design 
the model parametera A, B and m could be considered as correlated random 
variables representing the scatter in the data as well as modelling error. 

Ensuing reliability analysis, however, would be complicated. It is possible 
to significantly simplify the analysis by assuming that A is the only random 
variable. It will account for all material variability and modelling error; 

B and m are treated as constants. Estimates of the statistics of A are the 
mean = A and standard deviation = s. If A is assumed normally distributed, 
then R is a lognormally distributed random variable with statistics given 
by [10], 


Median of R = ^ = (A + B-p“) 

The coefficient of variation (COV) of R, 


C 


R 



/.434 


-1 


( 10 ) 


( 11 ) 


11 


Note that the master curve of Fig. 2 is in fact R. Also note, as stated 

above, that C = 12.3%, given on the figure, is a measure of the vertical 
R 

scatter. 

This form of computing ignores the statistical scatter of the esti- 
mates. However, it has been found that in the elementary least equares 
case, the distribution of the estimators does not play a significant role 
for sample sizes of roughly n > 30 [10], For the many available creep- 
rupture data sets of this size or larger, the statistical distribution 
of the estimates is not expected to be significant. 


BIAS AND UNCERTAINTY IN PREDICTING LONG TERM BEHAVIOR FROM SHORT TERM DATA 
The ability of a TTP to predict long time behavior (e.g., 20 to 40 
years) from short time data (e.g., less than 1 year) is of special interest 
to designers. A parameter which provides accurate predictions of long 
time properties would improve current design practices by reducing the 
risk and/or avoiding overdesign. The ability of a TTP to predict long 
time behavior can be measured by (a) constructing a model using short 
term data only, and (b) comparing long term data with the model predictions. 
Define A as the ratio of observed to predicted behavior; A is an index 
which defines the quality of a model. Following is a mathematical predic- 
tion of A in a form which is suitable for probabilistic design purposes. 
Define, 


R 


OBS(i) ” Observed strength (long term data) of i^^ specimen 

.th 


Rp„„,.x = Predicted strength of i specimen (model based on short 
term data only) 


n 


= niimber of long term data points 
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til 

Define A associated with the i specimen as. 


^ ^OBS(i) , i = l,n 
^ ®PRE(i) 


( 12 ) 


But A^, In general, will be different for each test specimen. Therefore, 
A will be a random variable. Compute the sample mean and sample standard 
deviation as. 


A 


In 
— - A 
n t X 

i=l 


(13) 




(14) 


It is assumed that A will have a lognormal distribution. In fact, statis- 
tical tests (not reported here) show that A seems to follow an approximate 
lognormal distribution for many data sets. The median and COV of a lognormal 
A will be respectively [10], 


'\j 

A 


= A//l + C2 


(15) 


C 


A ■ 



(16) 


Note that Cy^ measures both scatter in material behavior and strength 
modelling error. The closer A is to unity and the smaller is s^ , the better 
the predictions. If A < 1, then the model has a tendency to overestimate 
long time strength. For A > 1 the strength is underestimated. 
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Tables 2 and 3 summarize an exercise to study the behavior of A from 
typical creep-rupture data. Models based on the short term data are sum- 
marized in Table 2 for various data sets. The bulk of the short term 
data is less than 1000 hours, with only a few points between 1000 and 
3500 hours. Then the mean and COV of the bias in predicting long term 
behavior for each data set is summarized in Table 3. Lives of most of 
the long term data were from 12 to 18 months. The short term and long 
term data used here was provided by Goldhoff (7). Note that the MH param- 
eter seems to perform much better than the LM parameter as evidenced by 
the fact that A is consistently closer to one. 

RELIABILITY ANALYSIS 

In general reliability analysis, one should consider all sources of 
uncertainty including data scatter as well as the bias and distribution of 
modelling error. Material data scatter is described by the random variable 
A in the creep rupture analysis. But the bias A contains both scatter 
in material behavior and modelling error. In addition, temperature T 
and stress S can be treated as random variables reflecting uncertainties 
in the environment as well as the analysis. The goal of reliability 
analysis is to synthesize statistical information to compute risk. 

For a practical reliability analysis we need an estimate of the 
”in service” strength which accounts for data scatter as well as predictive 
capabilities of the model. Define the actual strength of the material 

R = f R 
a 


( 17 ) 



TABLE 2, Summary of the Analysis of Short Term Data 




































TABLE 3. Statistics of A 


MATERIAL 

LM I 


A 

C,% 

A 

INCOLOY 625 

0.477 

15.6 

UDIMET 500 

1 

_ . . 1 

0.928 

5.27 

UDIMET 500 

0.904 

7.69 

ASTROLOY 

0.884 

9.52 

RENE 41 

0.971 

6.22 

CR-MO 

0.967 

4.99 

HASTELLOY-X 

0.920 

15.6 

316-SS 

1.037 

8.25 

L-605 

0.668 

15.0 





Al-llOO 


1.036 


6.07 
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Lted by the TTP's. 


MH 

Number of Data Points 

A 

C,% 

A 

Short 

Time 

Long 

Time 

0.941 

12.00 

78 

17 

0.972 

3.95 

93 

12 

0.939 

6.93 

66 

37 

0.928 

8.65 

21 

12 

1.022 

8.19 

26 

11 

1.059 

8.98 

23 

10 

0.993 

19.0 

29 

18 

1.125 

9.45 

28 

10 

1.130 

17.8 

76 

28 

1.055 

7.38 

53 

11 
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where 't' is a random variable which accounts for bias and uncertainty of 
modelling error associated with using a TIP to extrapolate. R is the predicted 
strength of the material as calculated from the exponential model, Eq. 4, as 

R = log"J (A + BP”) (18) 


where P = P(t,T). R is a random variable because A and T are random vari- 
ables. Note that if T is assumed to be constant and A is normally dis- 
tributed, then R would have an exact lognormal. 

The statistics of ^ are established as follows. First it is asstamed 
that T = constant in the data used to compute A and . Then assume that 
¥ is lognormal. will have exact lognormal with statistics 

f\j ^ 

Median, R = R (19) 

COV, = /(I + c2) (1 + cj) - 1 (20) 


Now A contains both data 
Also note that 


scatter and modelling error so that C 


f\j *v 

R = A R (21) 

a 


Comparing Eqs. 19 and 21 and solving for in Eq. 20, the statistics 
for A are 


= 


f\j '\j 

Y = A 


A + 
/ 1 + 


- 1 


( 22 ) 


What we have done is to extract uncertainty due to modelling error from 
A . Thus , describes modelling error in the use of the TTP . 
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Consider a simple example which illustrates how to separate material 

variability from modelling error associated with using a TTP to extrapolate 

to long times. Use of the LM parameter on Incoloy 625 is demonstrated. 

is established from short term data as 0.0968 (See Table 2). From long 

term data (Table 3) , A - 0.477 and C. ® 0.156. The median of A is computed 

A 

*\j r\j 

from Eq. 15 as A = 0.471. Thus, ^ = 0.471. From Eq. 22, = 0.121. 

Thus, statistics on modelling error and material behavior are separated. 

It should be noted that C. > with equality when there is no modelling 
error due to extrapolation. The fact that exceeds in some of the 
data sets of Table 2 and 3, suggests that sample sizes were inadequate. 
Furthermore, note that contains some modelling error. The use of a 

TTP to describe a complex material phenomena suggests only an approxima- 
tion to physical reality. However, because the values of C (less than 10%) 

R 

do not significantly differ from those of other static properties, it is 

likely that this component of modelling error is small. 

If both the distributions of stress, S and R are lognormal, a closed 

a 

form expression for probability of failure is available [10] . The proba- 
bility of failure is 

p^ = $ (-S) (23) 

where $ is the standard normal cdf, and 3 is the safety index, as given by 

In (R /S) 

3 ^ 

/inia + c2)(l + c2^) 

S and Cg are the median and COV of S respectively. 


( 24 ) 
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In the more general case where T is also considered to be a random vari- 
able, reliability analysis becomes much more difficult- Techniques such as 
Monte Carlo or Rackwitz-Fiessler method [10, 12] must be employed to relate 
risk to the design parameters. 

It is important to note that computed values of p^ should be considered 
as ’’notional** values of the probability of failure- Because of the many 
uncertainties and assumptions and because risk levels are so low, it may be 
imprudent to argue that p^ defines risk levels in an actuarial sense. 

DEMONSTRATION OF RELIABILITY ANALYSIS 

Example 1. It is required to design a tension element to the following 
specifications (a) the temperature is 1100°F (593®C) , (b) the life of the 
component is 40 years (350,000 hrs) , (c) the applied load is lognormally 
distributed with statistics P = 10 kips (44.48 KN) and Cp = 25%. The materi- 
al considered is Hastelloy X. The problem is to find the component’s minimum 

4 

cross-sectional area A if the maximum allowable risk is p^ = 10 . From normal 
tables the target safety index is 3 = -$ ^(p^) = 3.72. In this problem 
both the IM and MH parameters will be used. The relationship of load to 

stress is S = P/A, hence ^5 ~ ” 25%. 

a) Larson-Miller Analysis 

The optimized form of the LM parameter for Hastelloy X is given 
by (See Table 2) 

P =(T + 46Q)(log^Qt + 18.59) 

In this example 

P = (1100 + 460) (log^Q 350,000 + 18.59) = 37.65 thousands 
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Using the data from Table 2 in Eq, 4, predicted strength is 

- 4.683 - 0.1082 (37.65)°*^^° 

“k = 25.48 ksi (175.7 MPa) 

The median and COV of the actual strength is obtained using the data 
of Table 3. Note that (a) both modelling error and material behavior 
are included and (b) the median of A is obtained using Eq, 15. 

= C. = 15-6% 

Ra A 

R =AR =0.909 • 25.48 ksi = 23.16 ksi (159.67 MPa) 
a 

Rearranging Eq. 24, ^ is given by 

R 

^ ^ ± 

exp {g[£,n(l + Cj^^)(l + Cg) ] ^} 

Then, 

S = 7.82 ksi (53.91 MPa) 

exp {3.72[£n(l + .156 ) (1 + .25 ) ] ^} 


For a safe design the minimum required area is 
A i I =^~2 = l*28in^(8.25 cm^) 


b) The Manson-Haferd Analysis 

The MH parameter for Hastelloy X is given by (see Table 2) 


T - 614 


87.8 


I log t - 11.08 

Using data from Table 2 in Eq. 4, predicted strength is 
logjj it * 2.905 - 3.777E-02 (87.8)®'®^^ 

R = 21.25 ksi (146.51 MPa) 
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The median and COV of the actual strength is obtained from Table 3. 

a. 

A is computed using Eq. 15. 

C„ = C = 19.0% 

Ra A 

= A R = 0.979 • 21.25 ksi = 20.8 ksi (143.40 MPa) 
a 

Then, as above, the maximum allowable stress is 

S = 5 — i— = 6.55 ksi (45.18 MPa) 

exp {3.72 [£n(l + .19 )(1+ .21'^)]'^} 

For a safe design, the requirement on the area is 

A i p/s = 10/6.55 =1.53 in^ (9.85 cm^) 

We are incapable of providing a commentary at this point regarding 
the poor agreement between the design as established by each parameter. 

From intuition, however, it seems reasonable to place more confidence in 
MH simply becuase it contains two (rather than one) empirical constants. 


Example 2. As an extension to Example 1, assume now that uncertainties 

in the temperature analysis are to be considered. It is estimated that 

the median (best estimate) of T is 1100®F and COV of T is 5%. The member 

2 2 

chosen has a cross sectional area of 1.60 in (10.3 cm ) , larger than the 
2 2 

1.53 in (9.85 cm ) required in the constant T case. Does the design 
satisfy the basic requirement that the target safety index, 8 > 3.72? 

A closed form solution is not available. To solve this problem, a 
numerical method, the Rackwitz-Fiessler (R-F) method will be employed 
[10, 12]. The computer program requires the limit state as input. 

(S = VR) = (S = 4' • log~J (A + B [(T + 460)(logj^Qt + 18.59)]“} 
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S, A and T are random variables. It Is clear that the computation 
of the probability of failure event would be most difficult by classical 
methods . 

The computations are summarized in Table 4. The computed safety 
index, 0 = 3.33, is less than the target safety index 0^ = 3.72, so the 
design is considered to be unsafe. It is interesting to note that if 
we let Cy 0 (T = 5^, a constant), then the safety index is 0 = 4.50. 

Clearly reliability is sensitive to uncertainties in T. 

We can also compute the notional probability of failure as 

Pf = $ (-3) 

For 0 = 3.33, = 4.36E-4, and for 0 = 4.50, = 3.44E-6. 

CONCLUSIONS 

A general method for reliability analysis of creep rupture data is 

presented. The method employs the time- temperature parameter (TTP) concept. 
An exponential model, log^^ R = A + B P™, relating strength R, to P, the 

TTP, was shown to provide a good fit to creep-rupture data. Evidence of 
the quality of the exponential model is provided by the low values of the 
coefficient of variation of R in Table 2 (typically less than 10%). Further- 
more, the exponential model fits well into a probabilistic design format. 

In general, the assumption that all of the statistical scatter of 
the data can be Itamped into A is valid only when the data is homoscedastic 
(constant scatter band) . This data is not always homoscedastic when the 
universal constants associated with each TTP are used. 
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TABLE 4. Summary of the Input and Output of the UA Reliability Program^^^ 
to Calculate 3, 


# Data 


Variable 

Distribution 

Mean 

Median 

cov% 

a 

S 

Lognormal 


6.25 

25.0 


A<W 

Normal 

4.683 


8.2 

0.382 

T 

Lognormal 


1100 

5.0 


f 

Lognormal 


0.909^*^^ 

13. 3^*^^ 



• Safety Index, B = 3.327 

• Notional Probability or Failure, = 4.40E-4 

• Design Point (approximately most probable value on the limit state function; 

close to the peak of the joint probability density function 
of the design factors) 

* 

S = 7.57 
A* = 4.674 
T* = 1329 
'F* = 0.857 


Notes 

(a) This program has the option of using the Hasof er-Lind, Rackwitz-^Fiessler 
or Chen-Lind method for computing the generalized safety index. 

(b) Hastelloy X; See Table 2. 

(c) ¥ = A; and A obtained using Eq. 15 and data of Table 2. 


(d) Computed using data of Tables 2 and 3 and Eq. 22. 
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The assumption that A is normally distributed is a matter of con- 
venience for the reliability analysis. While no formal goodness of fit 
test was used, the normal assumption for A seems to be very reasonable 
on the basis of visual inspection of R-P plots. 

Investigation of the bias in predicting long term behavior from 
short term data of both the LM and MH parameters suggest that both do a 
fair job. But no general trend was observed. Probability plots of A, 
and hence suggest that A (and follows a distribution close to a log- 
normal. 

Finally, failure is not always defined as rupture. It can be defined 
as excessive deformation. When this is the case, a similar approach to 
the one presented above is possible. Time to failure would indicate time 
to achieve a predetermined deformation, say 1, 2, or 5% creep. Creep 
strength would now be defined as stress, e.g., for 1, 2, or 5% creep 
strain. Fig. 4 illustrates the master curves for creep rupture and 1% strain. 



Fig. 4. Analysis of ^% Creep Strain Data. Creep Strain Data was Plotted by 


Using Optimized L-M Constant for Creep Rupture. 

flSTROLOY 


W 
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